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Abstract. In this paper we study a stochastic version of an inviscid shell 
model of turbulence with multiplicative noise. The deterministic counterpart 
of this model is quite general and includes inviscid GOY and Sabra shell models 
of turbulence. We prove global weak existence and uniqueness of solutions for 
any finite energy initial condition. Moreover energy dissipation of the system 
is proved in spite of its formal energy conservation. 



Introduction 

In recent years shell models of turbulence have attracted interest for their ability 
to capture some of the statistical properties and features of three-dimensional tur- 
bulence, while presenting a structure much simpler than Navier-Stokes and Euler 
equations. 

The main idea behind shell models is to summarize in a unique variable u n 
(usually complex- valued) all the modes with wave number k inside the shell A™ < 
|fc| < A™ +1 . Just like Navier-Stokes equations written in Fourier coordinates, the 
functions {m„}„gn satisfy an infinite system of coupled ordinary equations, where 
the non-linear term is quadratic and formally preserves energy. 

Shell models are however drastic modifications of Navier-Stokes equations. Firstly 
the variables {u„}„ e N representing three-dimensional shells (logarithmically equi- 
spaced) are one-dimensionally indexed by N. Secondly the shells are allowed to 
interact only locally. The choice to allow only finite-range interactions is a crucial 
simplification both from analytical and numerical perspective but it is well justified 
inside the Kolmogorov theory of homogeneous turbulence, where one neglects en- 
ergy exchanges between modes whose wave numbers differ for more than one order 
of magnitude. 

These two characteristics of the shell models represent both their weakness and 
their strength. The main weaknesses are the loss of the geometry and the restriction 
to questions concerning the turbulence energy cascade only. The strengths are sev- 
eral: from a numerical perspective, the lower number of degrees of freedom allows 
for more accurate simulations at high Reynolds numbers (although the implemen- 
tation of these simulations is not easy) ; from an analytical perspective, the simpler 
structure of the problem leads to sharper results both for the well-posedness and 
the understanding of the anomalous scaling exponents. 

A review of the subject that focuses in particular on these aspects can be found 
in Biferale [12] which is devoted to the turbulence energy cascade and collects 
results concerning the structure function S p (k n ) — E[\u n \ p ] together with numerical 
evidence and analytical conjectures about anomalous exponents. 

Main shell models. There are several different shell models in literature. The 
most studied are the GOY, introduced in Gledzer [22] and Ohkitani and Ya- 
mada [ ] and the Sabra introduced in L'vov et al. [ ]. Then there are two 
models with interactions that are somewhat simpler to study: one was introduced 
in Obukhov [28] and the other in Desnianskii and Novikov [20] and in Katz and 
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Pavlovic [23] . While all of the previous have variables indexed by N, there are also 
generalizations where the set of indexes is a regular tree (e.g. again in the first part 
of Katz and Pavlovic [23] and in Barbato et al. [3]) which is closer to a true wavelet 
formulation of Navicr- Stokes. 

For the viscous versions of GOY and Sabra, well-posedness, global regularity 
of solutions and smooth dependence on the initial data are known (Constantin, 
Levant and Titi [17]). On the other hand, for the inviscid case less is known, the 
state of the art being Constantin, Levant and Titi [18] where the authors prove 
global existence of weak solutions and, for sufficiently smooth initial conditions, 
uniqueness and regularity for small times. 

For the simpler shell models there are stronger results in the viscous case (see 
among the others Barbato, Morandin and Romito [9] and Cheskidov and Friedlan- 
der [13]), and even the inviscid case is understood quite well (the main results can 
be found in Katz and Pavlovic [23], Kiselev and Zlatos [24], Cheskidov, Friedlander 
and Pavlovic [14] and [15], Barbato, Flandoli and Morandin [5] and [8] and Barbato 
and Morandin [4]). 

Recently some stochastic shell models have been also proposed. An additive- 
noise version of the viscous GOY which is globally well-posed was introduced in 
Barbato et al. [ !]. The existence of invariant measures was proved in Bessaih and 
Ferrario [10]. In Manna, Sritharan and Sundar [27] and Chueshov and Millet [16] a 
small multiplicative noise version of the GOY model is studied; well-posedness and 
a large deviation principle are established. 

Finally, in Barbato, Flandoli and Morandin [(>] and [7] a stochastic version of the 
inviscid Novikov model was proposed, which is then generalized to the tree-indexed 
Novikov model in Bianchi [11]. In these last models the noise term is multiplicative, 
and it is tailored to be formally energy-preserving. The cited papers prove global 
well-posedness of weak solutions for both models and anomalous dissipation for 
the former. (By anomalous dissipation we denote the property by which the total 
energy of the system decreases in spite of the formal conservativity of the dynamics.) 

This type of noise is both elegant from an analytical point of view and physically 
meaningful in the sense that the interactions of Euler equations neglected in the 
shell models can be thought to be some sort of residual term which would behave 
(statistically) in a similar way. 

Main results of the paper. In this paper we study a general stochastic inviscid 
shell model, with a multiplicative noise term similar to the one in Barbato, Flandoli 
and Morandin [6]. We restrict ourselves to indexes in N, but we allow the variables 
X n to be multidimensional. 

This very general system of equations, given in (1), includes a stochastic ver- 
sion (51) of the inviscid GOY model (49) and a stochastic version (52) of the inviscid 
Sabra model (50). 

The noise term is formally conservative, in the sense that it acts only on the 
transport of energy, without giving or taking any part of it. One general way 
to obtain such a noise is the following. Suppose that some shell model satisfies 
j^u = B(u, u) with Re(u, B(u, v)) — 0, then also the system du — B{u, udt+aodW) 
will be formally conservative, the term a o dW being a perturbation acting only on 
the transport. 

The first aim of this paper is to prove that for the general model (1) there are 
global weak existence and uniqueness in law of I 2 solutions. This result, to the 
knowledge of the authors, is the first result of global well-posedness for the inviscid 
GOY and Sabra models, both deterministic and stochastic. 

The second important result concerns anomalous dissipation. Theorems 7 and 
8 state that for both stochastic inviscid GOY and Sabra models energy decreases 
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with positive probability at all times, that it becomes arbitrarily small again with 
positive probability, and that if the initial energy is small enough, the solution 
converges to zero at least exponentially fast almost surely and in L 2 . 

Strategy and organization of the paper. Section 1 introduces and describes 
the general stochastic inviscid shell model (1). This is a subtle matter, since the 
requirement that the noise acts on the transport term in a conservative way leads 
to several algebraic conditions. 

One of the key ideas of the paper is to use Girsanov theorem to study the 
original problem through an auxiliary linear system. Sections 2 and 3 are devoted 
to establish the relation between the linear and non-linear systems. 

In Section 4 from the linear system we deduce the evolution equations for the 
second moment of components, which turn out to form a deterministic linear system 
and can be conveniently studied through the theory of q-matrices of continuous- 
time Markov chains. The first consequence is uniqueness of solutions: strong for the 
auxiliary linear system (Theorem 2), and in law for the main model (Theorem 3). 

Existence of global solutions is classical and straightforward, and is detailed in 
Section 5. This concludes well-posedness for the main model. 

Section C is devoted to anomalous dissipation, which is deduced from the be- 
haviour of the continuous-time Markov chain associated to the equations for the 
second moment of components. In particular the chain turns out to be dishonest, 
in the sense that a.s. it reaches infinity in finite time. This "loss of mass" pulled 
back to the initial system becomes a loss of energy towards higher and higher com- 
ponents, which is what we call anomalous dissipation. 

To formalize the link between the chain and the main model, one needs two steps: 
Borel-Cantelli lemma to get a.s. statements about the auxiliary linear system, and 
Novikov condition for Girsanov theorem to pass from the auxiliary system to the 
main model. Theorem 6 shows that if the initial energy is small with respect to the 
noise, Novikov condition holds also at t = oo, and so even the exponential decay of 
energy can be deduced. 

Finally we need to show that indeed GOY and Sabra models can be included 
in the general model (1) and summarize the results for these two models. This is 
done straightforwardly in Section 7. 

1. Main model and formal requirements 

The general model of this paper is equation (1) below. Since it is both complex 
and written in a synthetic but unfamiliar way, it will be helpful to start with a 
particular example to be kept in mind as a reference. 

Consider on a complete filtered probability space (f2, J : 00l {jFt)t>o, P) the infinite 
system of stochastic differential equations in Stratonovich form below 

dX n — \ n ^ X 2 x _^dt — \ n X n X n +idt 

+ X n ~ 1 X n -i odW„_i - \ n X n+1 odW n , n > 1, 

where (W n ) n >o is a sequence of independent Brownian motions and Xq = 0. This 
model is a stochastic version of the inviscid Novikov shell model, and was introduced 
in Barbato, Flandoli and Morandin [6] and [7]. The two deterministic terms are 
coupled in such a way to cancel when we sum X n dX n over n. Apart form this, they 
are of the same form, representing an interaction between X n and the product of 
two other components. The two stochastic terms are coupled among themselves in 
the same way, and moreover each of them is associated to one of the deterministic 
terms, so that the equation rewrites 

dX n = X 1 ^ 1 X n -i(X n -idt + odW n ^) - X n X n+1 (X n dt + odW n ), n > 1. 
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We want to generalize this model to different types of interactions and multidi- 
mensional structure. We will also try to keep the notation as less cumbersome as 
possible, and in this view we will rewrite this as a sum over a set of interaction 
terms /, that will include pairs of cancelling interactions. 
We are finally able to write the general model. 

(1) dX n = h,nBj (X n+ri , X n+hi dt + a orfWj, n +; t< ) , n > 1. 

iei 

Here each X n is a d-dimesional real- valued stochastic process, I is some finite set 
with an even number of elements and (Wi in )ig/.„ e z is a family of ci-dimcnsional 
brownian motions (independent apart from some deterministic relations explained 
below). For all n e Z and i E I, ki, n is a real constant, Bi is a bilinear operator on 
R d while ri and hi are integer numbers. 

In the example of Novikov model given above, I has two elements 1 and 2, d = 1, 
Bi(a, b) = ab for i = 1,2, the coefficients are given by 



i 




1 


-1 -1 A"- 1 


2 


1 -A" 



and the Brownian motions are independent apart from W\^ n = W2, n a.s. for all n. 

Going back to the general model, since and hi may be negative, we pose 
X n = for n < and = for i, n such that n + ri < or n + hi < 0. We will 
also require that h := maxi hi > otherwise X n is constant for n < —h. 

We now list a first set of basic requirements on these models 

i. Finite range: I is a finite set. 
ii. No self interactions: ri ^ for all i € I. 

in. Exponential coefficients: ki, n = X n ki for all i G /„ and n > 1; here A > 1 
and fcj are real numbers and /„ := {i £ I : n + ri > l,n + hi > 1}. If 
i £ I n then fcj jn = 0. 

The fourth but very important requirement is the formal (also called local) conser- 
vation of energy, which is assured by some cancellations, as described below. The 
intuitive meaning of the conditions detailed below is that / must be formed by pairs 
{i, i} of cancelling interactions such that for all n there exists h = n + ri such that 

h,n (X n , Bi(X n+ri ,X n+hi dt + a odWi in+ /,.)) 

+ ki <n (X n , Bi(X, h+ri ,X, h+hl dt + a odW^ n+h -,)) = 0. 

Here (•, •) denotes the scalar product in R d . 

To make things formal and clean we need to introduce this definition. 

Definition 1. Suppose r is a permutation of / with no fixed point and such that 
t = r _1 . Let /* be any subset of / such that / is the disjoint union of /* and 
t(I*). We say that a family of d-dimensional Brownian motions W = (Wj in )j e j ine z 
is symmetric with respect to r if the restriction of W to I* x Z is a family of 
independent Brownian motions and W T ^ in — W^ n a.s. for all i £ I and n€Z. 

Clearly this definition does not depend on the particular choice of I* , but nev- 
ertheless by invoking this definition, we will implicitly suppose that we are fixing 
the set I*. 

We can now state the fourth requirement. We will use the notation i = t(i). 



STOCHASTIC INVISCID SHELL MODELS 



5 



iv. Local conservativity: there exists r such that W is symmetric with respect 
to t and the following relations hold for alH £ / 

(2) fcj = -fc i A" ri , 

(3) (u,Bi(v,w)) = (v,Bi(u,w)), Vti,D,mel' ! , 

(4) n = -n, 

(5) hi = hi-ri, 

Remark 1. These algebraic requirements are meaningful, in the sense that given 
I and t, there exist k, B, r and h satisfying them. Truly, it is easy to verify 
that however we define these objects on /*, there is exactly one extension on all I 
satisfying the above conditions. 

The following lemma summarizes some other trivial but useful consequences of 
the requirements above. 

Lemma 2. Let ip be the automorphism on I x Z defined by ip(i,n) := (l,n) :— 
(r(i),n + Ti). Then there exists A C I x Z such that ip(A) = A c . Moreover the 
following relations hold for all i £ I and n G Z 

k- - — — fr- 
n = n + r.; 
n = n + rj 
W^n +hi = W^ n+hl a.s. 
In particular it is now straightforward that for all i and n, 

(6) ki, n (X n , Bi(X n+ri , X n+hi dt + a odW^n+hi)) 

+ h.n {X? l ,Bi{Xf l+r% ,Xf l+ h-,dt + a odWifi+h*)) = 

Since by the Stratonovich form of Ito formula we have 

d(X n ,X n ) = 2{X n , odX n ), 

if we sum formally these quantities over n substituting (1) and using (6), we have 
J2 n >i d\X n \ 2 — 0, so we may expect these models to be conservative. 

Actually, this is in general not true. Rigorous arguments in the following sections 
will show that dJ2 n >i \X n \ 2 ^ J2 n >id\X n \ 2 anc ^ that J2 n >i \X n \ 2 decreases with 
positive probability. 

2. Ito formulation and auxiliary equation 

We prefer to reformulate equation (1) with Ito integration. Proposition 3 below 
states that the equivalent Ito differential equations are the following 
(7) 

dX n = ^ k l ,7 l B l (X n+ri , X n+hi dt + adWi, n+ hi) ~ y X! k ln L i x ndt, n>l 
iei iei 

where Li is the linear map on M d given by Li :— BiBf . (Here Bi is interpreted as 
a linear map from R d to R d .) In components, L"' 13 = s Bf n ' A 'B^' 1 ' 6 . 

We will also introduce the auxiliary linear system of equations and their solu- 
tions. This will be needed afterwards, for Girsanov theorem. 

2 

(8) dX n = h,nBi{X n+n ,adWi, n+hi ) - — ^ k 2 n LiX n dt, n > 1 
iei iei 

Let H := Z 2 (M d ) denote the state space and || • || its norm. 
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Definition 2. Given an initial condition x £ H, a weak solution of non-linear 
system (7) (respectively of linear system (8)) in H is a filtered probability space 
(fi, Too, (-?t)t>0) P)i along with a family of Brownian motions W, and a stochastic 
process X such that 

i. W = (Wj lTl )i£/ in gz, is a family of <i-dimensional Brownian motions on 
(fi,.^, (J r t )j>o,i :> ) adapted to the filtration and symmetric with respect 
to t; 

m. X = (X„)„>i is an _ff-valued stochastic process on (f2, .Foe (J-t)t>o, P) 

with continuous adapted components; 
ni the following integral form of non-linear equation (7) holds for all n > 1 
and all t > 0, 



(respectively, the following integral form of linear equation (8) holds for 
all n > 1 and all t > 0) 



The next proposition shows that what is defined above is actually a solution of 
the Stratonovich formulation of the non-linear system. 

Proposition 3. If X is a weak solution of the non-linear system (7), the Stratonovich 
integrals 



(9) 






kl n L t X n (s)ds} 



(11) 




are well defined and equal to 



(12) 




Hence X satisfies the Stratonovich equations (1). 



Proof. We write components explicitly, in particular Bi = (_B[ 
s Bf' 1 ' 6 B^ n ' S . Component a of (11) rewrites 



and L° 



(13) 




The stochastic integral can be rewritten in Ito form 



(14) 




(s) odWl 



i,n-\-h. 




i,n-\-hi 



J f 
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we only need to compute the quadratic covariation term. When n + < 0, this is 
zero, while if n + ri > 0, by writing (9) with n + ri in place of n, we find 

(15) [x? +ri ,W? n+hi ] t 

E E B i / ^i+ri ^M-n+fc, 



jei 6, v Ja 

The very last quadratic covariation differential can be (is or 0, depending on whether 
the two particular BM involved are equal or independent. They are clearly inde- 
pendent when j (jL They are independent also when j = i, since r% 7^ 0. 
Finally, they are a.s. equal when j = % and 77 = 7 by conditions (5) and (6). We get 



(16) 



X 



w, 



7 



n+ri > " v i,n+hi 



S Jo 

= -ah, n J2Bt^ fx s n {s)ds 
x Jo 



where we used conditions (2), (3), and (4). Putting all together 
(17) 

Tk itn Bi(X n+n (s), odW lyn+hz (s)) - / <7k itn B i (X n+ri (s) 1 dW l , n+hl (s)) 



0,1 

t „2 



X 



w 1 

n+ri ' i,n+h 



tCE^ 7 E 5 



<5,/3,7 



0,1 



6 

t ~2 



Jo * s p n L/o * 

The latter is correct also when n + r"j < 0, since in that case fcj n = 0. 



X S n (s)ds 



(s)ds 



□ 



Definition 3. Given an initial condition x 6 H, an energy controlled solution of 
the non-linear system (7) or the linear system (8) is a weak solution of the same 
system of equations in the class L°°(f2 x [0, 00); H). In particular, if = ||a;||, 

it is called a Leray solution. 



3. GlRSANOV TRANSFORMATION 

We turn our attention to the terms X n+ h.ds + adW^n+hi m equation (7). We 
would like that, under a new probability measure Q, these were the differentials 
adYi^ n+ hi , where Y is again a family of d-dimensional Brownian motions symmetric 
with respect to r. To do so, we state an infinite-dimensional version of Girsanov 
theorem whose proof can be found in Da Prato et al. [19]. 

Theorem 1. On a filtered space (Q, Too, (Ft)t>n,P), let (Wj)jeN be a sequence of 
1- dimensional adapted independent Brownian motions and let X = (Xj)j e iq be a 
sequence of adapted semimartingales such that ■ Xj(t) < 00 for all t > 0. 

Let Yj(t) := /„' Xj(s)ds + Wj(t) for j € N. Put, for < t < 00 

Mt = expj- jf £ X, ;( S )rfW, \a) - \ jf E *i («)<**} 
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Fix < T < oo. Suppose E[Afy] = 1, i/iera M zs a closed martingale on [0, T] 
and ffte density ^ — Mt defines a new probability measure Q on Tt under which 
(Yj-)jgN is a sequence of independent Brownian motions on [0, T). 

Moreover, to prove that E[M<r] = 1, Novikov condition can be used, namely it is 
enough to prove that 



(18) E 



exp{^ T ][>J( S K< 



< oo 



By virtue of Theorem 1 it is quite easy to verify that, by changing the prob- 
ability measure and the family of Brownian motions, any Leray solution of the 
non-linear system (7) can be transformed into a Leray solution of the auxiliary 
linear system (8). The reverse is also true. 

In our notation, X will always denote the solution process; P and Q will denote 
the probability measures for the non-linear and linear systems respectively; Pt and 
Qt their restrictions to Tt\ W and Y will denote the two associated families of 
Brownian motions. 

The relation between W and Y is ensured by the following definition. For i 6 J, 
n G N and t > 0, let 

r* i 

(19) y i>n (i = / -X n {s )ds + W i>n (t) 

Jo & 

Suppose we can define the two martingales 

(20) Z t = - f Ysi^^+hM^dW^+hM), 

(21) Z f = f y2(a- 1 X n+hi (a),dY i ,n+h i {8)) 

Then it easy to verify that 

Z t = -Z t + \ ( V \X n+hi (s)\ 2 ds = -Z t + [Z,Z} t = -Z t + [Z,Z} t 

nGN 

so that Z t - \[Z, Z] t = -Z t + \[Z, Z] t . We will then pose 

(22) ^ = exp{Z t - \[Z, Z] t ], ^ - exp{Z t - ±[Z, Z] t } 
We are now ready to make a precise statement 

Proposition 4. Suppose (f2, J^oo, (Jtjoo, P, W, X) is a Leray solution of the non- 
linear system (7). Fix any < T < oo. Let Qt be the measure on Tt defined 
by (22) and let Y be the family of Brownian motions defined by (19). 

Then Qt is a probability measure and (n,J-T,(J~t)o<t<TiQT,Y,X) is a Leray 
solution of the linear system (8) on [0, T] . 

Proof. By (19), equation (9) is equivalent to equation (10) with W replaced by 
Y. We need to prove that Qt is a probability measure and that Y is a family of 
Brownian motions symmetric w.r.t. r, hence we apply Theorem 1. 

The sequences we use are —X^ +h . and W- n+h ., both for 1 < j < d, n G N and 
i G I*. (Notice that we use I* instead of I since we need the independence of the 
Brownian motions.) 
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By Leray property and the finiteness of I* , we get a very strong bound 
(23) J2J2\X n+ht \ 2 (t)<\r\\\x\\ 2 a.s. for all t > 



iel' n>l 



by which we immediately deduce both that Z t in (20) is well-defined and, by the 
finiteness of T, that 



[Z, Z] T = [Z, Z] T = f T E \\Xn +hl \ 2 (s)ds < ™ 



x\\ 2 T 

2 ' 



(24) 

hence Novikov condition holds, namely 

(25) E[e^ z H = E[e^H < exp \I3^L < oo 
Finally, Y is symmetric w.r.t. r since both P-a.s. and Qt~&.s. 

(26) Y i>n (t) = f -X n {s)ds + W t , n (t) = f -X n {s)ds + W itn (t) = Y^ n (t) □ 

Jo c Jo ° 

The converse is also true. We give it without proof since it is almost identical to 
the one above. 

Proposition 5. Suppose (J7,^"oo, (Ft)t>0iQ,Y,X) is a Leray solution of the linear 
system (8). Fix any < T < oo. Let Pt be the measure on Tt defined by the 
second one of (22) and let W be the family of Brownian motions defined by (19). 

Then Pt is a probability measure and (n,J-T,(J~t)o<t<T,PT-,W,X) is a Leray 
solution of the non-linear system (7) on [0,T]. 

Remark 6. By Caratheodory theorem, the family of probability measures (Qt)t>o 
extends in a unique way to some probability measure Q on and the same stands 
for (Pt)t>o- Hence solutions of non-linear and linear systems are associated also 
on infinite time span. From now on we will drop the t and use the symbols P 
and Q with this meaning. One should anyway keep in mind that while Pt and Qt 
(and hence P and Q) are equivalent on Tt for any finite T, they are not in general 
equivalent on J-^. 



4. Closed equation for E q [|A"„(i)| 2 ] and uniqueness 

Denote by the mathematical expectation on (f2, J 7 , Q). It turns out that if 
Li is the identity for all i £ I, then E^ [|X n (£)| 2 ] satisfies a closed linear differential 
equation which will shed new light on the behaviour of solutions, in particular by 
giving an easy way to prove uniqueness. 

Proposition 7. Suppose (f2, .Foo, (J c t)t>o I Q-, Y, X) is an energy controlled solution 
of the linear system (8) . 

Then for all n > 1 and t > 

(27) jE Q [\X n \ 2 ] = -Y,^kt n E Q {X n ,L^ n )+J2^kt r fi Q (X n+rt ,Lr l X n+rz ) 
iei ieJ 

Proof. We start by computing the quadratic variation of X n . We use (10) and 
the independence of Y^n+foj and Y^ n _|_/j. when j =/= i. (If j = i, then n + hj = 
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n + hi — T\ 7^ n + hi.) 
d[X n , X n ] t 

iei je 



^^h^k^nd / Bi(X n+n ,dY^ n+hi ), I B 
-iei iei L^o Jo 



— 'y ' fcj n <^ / Bi{X n j rTi , dYi jn+ hi)> / Bi(X n+ri ,dYi^ n+ hi 
iei Uo Jo 

Eh 2 R Q,/9:7 R a:<5,e v/3 jrc5 j [yy y£ 1 

= X] ^> 5f' a,1 Bf t>,7 lf +r .I* +r ii 



3 (^n+rj ) ^^j',n+/ij ) 
) 



^i,n (-^n+r< 7 LjX n+ri )dt 



iei 



We also used (3) and the definition of = ^ 7 s B"' 1 ' B 
Now we are able to compute the differential of |X„(t)| 2 

(28) d(\X n {t)\ 2 ) = 2{X n (t),dX n (t))+d[X n ,X n ] t 

= 2^ak hn (X n ,Bi(X n+rt ,dY itn+hi )) -J2< j2k tn( x n,L i X n )dt+ 



iei 



>ei 



iei 

By the definition of energy controlled solution, |X„(i)| < ||X(t)|| < C almost surely, 
so m particular E Q f* \X n \ 2 (s)\X n+ri \ 2 (s)ds < oo and hence the first term above is 
a martingale, with mathematical expectation zero. If we now take the mathematical 
expectation of the integral form of (28) we get the integral form of (27) and we are 
finished. □ 

When all the Li are the identity, equation (27) becomes a linear closed differential 
equation 

(29) -V [\X n \ 2 ] = -XVfc? n E« [|X„| 2 ] +J2^tn^ Q [\Xn + rf] 



iei 



iei 



We notice that this system of equations is of a peculiar kind, with negative di- 
agonal and non-negative off-diagonal entries, thus suggesting a connection to the 
Kolmogorov equations for continuous-time Markov chains on the positive integers. 
We investigate this relation presently. 

Denote by II = (7r m , n ) m , n >i the infinite matrix associated to this system: for 
n, m > 1 and m ^ n let 



(30) 



'* n , m := £ a 2 k 2 n = a 2 X^ £ k 2 
iei iei n 

ri=m—n ri—m—n 

n,n '■= ~ £ <j2 ^'i,n — — cr 2 A 2 ™ £ kf =: 



iei iei n 
Remark 8. Recall that for n > 1, I n := {i € / : n + ti > 1, n + hi > 1}, so that 

I\ C I2 c • • • C I no = I no +l = ■■•=/. 
where no > 1 — min{ri, hi : i £ I}. Hence for example n n = 0(X 2n ) as n — > 00. 
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Corollary 9. Suppose that (Q,X) is an energy controlled solution of equation (8) 
and that all the Li are the identity. Then u — (u„)„>i defined by u n (t) — [|X„(i)| 2 ] 
is a non-negative solution of the Cauchy problem 



v! = uli 

u n (0) = \x n \ 2 n > 1 



(31) 

in the class L°°([0, oo); Z 1 ) . 

Proposition 10. The infinite matrix II defined above is the stable, conservative 
q-matrix of a continuous-time Markov chain on the positive integers. Moreover H 
is symmetric. 

Proof. II is a stable q-matrix if 7r n .„ < for all n, 7r„. m > whenever n ^ m and 
for all n 

(32) 2_j 7r ™ 1 " 1 - 7r ™' 

moreover it is conservative if the latter holds with equality. 

The first two conditions are obvious, and the third one follows from the fact 
that the sets {i € I : n + r< = m} with m > 1, m ^ n form a partition of 
{i e I : n + n > l,fei,„ ^ 0}. 

Finally, for m =/= n, 

(33) TT n . m = ^ °" 2 *i,n = °" 2 **>-n = V 2k lm = Tm,n □ 



iei iei iei 

—m—n n=n—m Vi—n—m 



The q-matrix of a continuous-time Markov chain is associated to the forward 
and backward Kolmogorov equations, namely 

(34) u' = ull 

(35) vl = Tlu 

The transition probabilities of the Markov chain p n , m (t) solve both equations, in 
the classes l l and l°° respectively, with fixed n and m respectively and with initial 
condition w„, m (0) = 5 n ,m- 

These equations always have at least one shared "special" solution fij{t), which 
is a transition function, and is called the minimal solution. They do not always 
have uniqueness of solutions. Here it will be important that there is uniqueness for 
the forward equation and not for the backward. The key information is that the 
q-matrix is symmetric. 

Lemma 11. Suppose II is a stable and symmetric q-matrix. Consider the for- 
ward equations with zero initial condition. Then the only non-negative solution in 
^([Ojoo);* 1 ) is zero. 

More in general, given any non-negative I 1 initial condition, there a unique so- 
lution in the same class. 

Proof. For the first part, we follow the classical approach by Laplace transform, 
introduced by Feller [21]. Let p be such a solution. For all n > 1 and t>0we have 

' P'rSf) = J2kPk( t W,n 
Pn(t) > 
Pn{0) = 

lE k Pk(t)<c 

For all n > 1, let z n = / °° e~ t p n (t)dt. Clearly J2 n Zn — ^> so we can choose m 
such that z m > Zk for all k. 
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Notice that, since II is stable and symmetric 



m 



(t) + TT m C < 2Civ m < oo 



hence we can integrate by parts and use symmetry and stability again to get 



We conclude that z n = and p n = for all n. 

For the general case, let fi,j(t) be the minimal solution of II and let u° be a 
non-negative, I 1 initial condition. Then u n (t) = J2i>i u ifi,n(t) is a solution in the 
required class. Let v be another such solution and let p = v — u. By a forward 
integral recursion (FIR) approach it is easy to show that the minimality of / passes 
to u, in that v n (t) > u n (t). (See for example Anderson [I], Theorem 2.2.2.) So p is 
a solution of the same problem, but with null initial condition and the first part of 
the lemma applies. □ 

We have finally collected all elements to prove uniqueness of solutions. 

Theorem 2. Suppose Li is the identity matrix for all i G /. Then there is strong 
uniqueness for the linear system (8) in the class of L°°(Q x [0, oo);H) solutions. 

Proof. By linearity of (8) it is enough to prove that when the initial condition is 
x = there is no non-trivial solution. Suppose (Q,X) is any energy controlled 
solution with zero initial condition, then by Corollary 9, Proposition 10 and the 
first part of Lemma 11, [|X„(i)| 2 ] = for all n and t, hence X — a.s. □ 

Remark 12. This result applies seamlessly also to the case of L°°, non-anticipative, 
random initial conditions. 

Uniqueness of solutions for the auxiliary linear system is then inherited by the 
original non-linear system, but in a weakened form. 

Theorem 3. Suppose Li is the identity matrix for all i € / and let T > 0. Then 
there is uniqueness in law for the non-linear system (7) in the class of Leray L°°(f2x 
[0,T];H) solutions. 

Proof. Suppose we are given two solutions (pW , W {1 \ X^) and (p( 2 \ , X^). 
We want to prove that 



(36) E p,1, [/(I< 1 )(i 1 ) I l( 1 )(t 2 ) 1 ..,l( 1 >((„))] 

= E pl2, [f(X^(t 1 ),X^(t 2 ),...,X^(t n ))} 



where / is any bounded measurable real function on H n and ti,t%, ■ ■ ■ ,t n G [0, T]. 
By Proposition 4 and the first one of (22) we have that, for j = 1, 2 



(37) E pU> [/(!»(*!), X^(t 2 ), . . .,X^(t n ))] 

= E« U) [exp{-4 j) + \\Z^\ Z«] T }/(X«(^),X«(t 2 ), . . .,X^(t n ))] 



Where Z^> is defined by (20). By Theorem 2 equation (8) has strong uniqueness, 
hence we can apply an infinite-dimensional version of Yamada-Watanabe theorem 
(see Revuz and Yor [31] or Prevot and Rockner [ ]) to deduce that the laws of 





k^rn 



k^m 
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(X^\W^) and (X^ 2 \W^) on C([0, T]\ M 2rf ) N are equal, under QW and 
respectively. 

Then of course we can include also Z^ by their definition and conclude that 
(X^,W W , Z (1) ) under QM and (X< 2 ), W (2 \ Z (2) ) under Q( 2 ) have the same law, 
yielding in particular that (37) does not depend on j. □ 

5. Existence of solutions 

In this section we prove strong existence of solutions for the linear auxiliary 
model and deduce weak existence for the non-linear model. The approach is by 
finite-dimensional approximation and follows Pardoux [29] and Krylov and Ro- 
zovskh [ ]. The border term of the finite-dimensional systems is chosen so that 
energy conservation holds, giving a strong tool to prove convergence. 

Theorem 4. Let J 7 ^, (J r t)t>o, Q) be a filtered probability space. Let Y be a 
family of adapted d- dimensional Brownian motions symmetric with respect to r. 
Given an initial condition x € L°°((Q,Fq);H) and T > 0, there exists at least 
an H-valued stochastic process X with continuous adapted components, such that 
Q-almost surely \\X(t)\\ < \\x\\ for all t > and for all n>l and all t>0, 
(38) 

X n (t) = Xn + zZ{j Q <7ki,nBi(X n+ri (s),dWi, n+hi (s)) - y fcf„i l X„(s)ds| 

Such a process is called strong Leray solution. 

Proof. For every positive N, let Am = {1, . . . , N} and consider the finite dimen- 
sional stochastic linear system 
(39) 

'dxi N) = J2 ieI k hn B l (X i n 1 l ) ri: adY hn+hi ) - ^ E*e/ kf^X^dt, n £ A N 
Xi N \0) = X n, neA N 
xi N) = Q, n€Z\A N 

This system has a unique global strong solution X^ N \ By the local conservativity, 
we can prove that the I 2 norm is Q-a.s. constant. In particular we notice that 
equation (28) applies to X^ N \ for n g An without modifications. Then we sum on 
n and apply Lemma 2 to get 

£ d(|lf)(t)|^2^< n (lf)B 1 (4!l^,)) 

n£A N JxZ 

- J2 c- 2 k 2 n (xi N ) lLl x^)dt + J2 «ia,i^))rft 

JxZ JxZ 

= 2aY,{k,n(X( N \B i (xi J ll i ,dY i>n+hi )) + ki fi (xi N \Bi(xi 1 ll v dY, fi+h A) 
- J2 ^kl n {Xi N \L iX W)dt + £ a 2 kUX^\ Li X^)dt = 

JxZ JxZ 

Thus 

(40) Yl \ X i N) (*)\ 2 = £ IXn| 2 <IMIi-(n ;H ), V*>0, Q-a.s. 

meaning in particular that the sequence X^ N > is bounded in L°°(f2 x [0, T];H). 
Hence there exists X in the same space and a sequence N/~ f oo such that X^ Nh ' 
X as k — > oo. A fortiori there is also weak convergence in L 2 (fi x [0, T]; H ). 
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Let X denote the subspace of L 2 (Q, x [0,T];H) of the progressively measurable 
processes, then G X for all N. The space X is complete, hence it is a closed 

subspace of L 2 in the strong topology, hence it is also closed in the weak topology, 
so X must be progressively measurable. 

Now we want to show that X indeed satisfies equation (10) with Y in place of 
W, as each of the I^'s docs. Fix t G [0,T], i G I, n > 1 and l,m G {1,2, ...,d}. 
The map 

[ Vl +ri {s)dY™ +hi {s) 
Jo 

is a linear strongly continuous map from X to L 2 (il; R), so it is weakly continuous. 
Equation (10), written component- wise, reduces to a finite sum of one-dimensional 
stochastic integrals like the one above, hence we can pass to the limit and so X 
solves the same equations. A posteriori, from these integral equations, it follows 
that there is a modification such that all components are continuous. 

Finally we prove the Leray property. Consider the product measure fi = Q x C 
on O x [0,T]. Let e > 0, let A = {(u),t) : \\X(u),t)\\ > ||x(w)|| + e} and let 
U = pfjflA € L 2 {tt x [0,T];iJ). Then 

(X,U) L *= J \\X\\l A dn>J \\x\\dn + en(A) 

On the other hand by Cauchy-Schwartz inequality on H and by (40). 

(X^\U) L2 = J(x^\^\ H l A dn< I \\X^\\l A dn< J \\x\\dfi 

Taking again the weak limit, we have fi(A) — and by the arbitrarity of e, we get 
that /i-a.e. ||X|| < ||x||. This can be improved by the continuity of components, 
which implies that the maps t >— > ^2 k<n \Xk(t)\ 2 are continuous and hence Q- 
a.s. bounded for all t and all n by ||xl| 2 - Letting n — > oo we conclude. □ 

The strong existence statement for the linear model becomes a weak existence 
statement for the non-linear model, due to Proposition 5. 

Corollary 13. Given an initial condition x G H and T > 0, there exists at least 
one Leray solution of the non-linear system (7) in the class L°°(rt x [0,T];H). 

6. Anomalous dissipation 

In this section we want to prove that ||X(t)|| goes to zero in some sense. We con- 
sider the differential equation for the second moments (29) and study the continuous- 
time Markov chain that has it as its forward Kolmogorov equation. The following 
proposition gives an explicit connection between the two. 

Proposition 14. Suppose Li is the identity matrix for all i £ I. Let x G H and let 
(<5, X) be the unique Leray solution of the linear system (8) with initial condition x. 
Then there exists a continuous-time Markov chain (£,t)t>o defined on a probability 
space (S,S,V) taking values in N and with q-matrix LI defined by (30), such that 
for allt>0 

(41) EQ[\X n (t)\ 2 ]=\\x\\ 2 r(^=n), Vn>l 

(42) ^ Q [\\X(t)\\ 2 ] = \\x\\*Pfo G N) = \\x\\ 2 V(r > t) 
where 

t := sup{£ : £ has finitely many jumps in [0, t)} G (0, oo] 
is the so-called explosion time of the Markov chain. 
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Proof. Let p° n := \x n \ 2 /\\x\\ 2 for all n > 1. It is standard to formally construct 
a continuous-time Markov chain £j on (S, S, V) with initial distribution p° and 
rates 7r„ jm as defined in (30). Heuristically, the process starts at a random position 
£o with V(£o = n) = Then every time the process arrives in a position n it 
waits for an exponentially distributed random time with rate 7r„ = ^2 m ^ n Kn,m and 
then jumps to a new random position different from n cheosen with probabilities 
( n "' m ) ^ ■ This defines £t up to time r. At time r we say that £ has reached the 

boundary. (Sometimes this is done by adding one absorbing point 9 to the state 
space.) 

For n > 1, t > 0, let p„(t) := 7-'(^t = Then p is a non-negative solution of 

iy = P n 

\p(0)=p° 

in L°°([0, oo); Z 1 ). By Corollary 9, u/||x|| 2 is another such solution, hence by the 
uniqueness result in Lemma 11 we have proved (41) and by summing up also (42). 

□ 

The following is the main result for the anomalous dissipation of the auxiliary 
linear system. The exponential decay of the expected value of energy follows from 
the Markov property of the chain. 

Theorem 5. Suppose Li is the identity matrix for all i 6 /. Let x G H and let 
(Q, X) be the unique Leray solution of the linear system (8) with initial condition 
x. Then the quantity [|| X(t) || 2 ] is strictly decreasing in t . Moreover there exists 
a constant /i > depending only on the coefficients ki, n and a constant C > ||x|| 2 
depending only on fcj, n 's and x, such that for all t > 

EQ[\\X(t)\\ 2 ] <Ce-^ 

Proof. By Proposition 14 we are given a probability space (S, S, V) and a continuous 
time Markov chain £ on the positive integers, defined up to some stopping time r 
such that TE,Q[\\X(t)\\ 2 } = \\x\\ 2 V(t > t), so we study the latter probability. 

Once we will prove the second statement, the fact that V(r > t) is strictly 
decreasing in t will follow from V(t = oo) < 1 by Chapman-Kolmogorv equation. 
(This is a standard exercise on continuous-time Markov chains whose proof is not 
difficult. See for example Lemma 13 in Barbato, Flandoli and Morandin [7].) 

We want to prove the exponential bound. Let £fc for k = 0,1,2,... be the 
discrete time Markov chain embedded in £, meaning that = £ t for t between the 
fc-th and the (k + l)-th times of jump of £. 

For n > 1, let V n := > 1 : Cfe = n} be the number of times £fc visits n. 
The law of V n , conditioned on V n ^ is geometrically distributed. Since the sum 
of a geometrically distributed number of i.i.d. exponential r.v.'s is exponentially 
distributed, the total time T n spent by £t on n, conditioned on ever reaching that 
site, is exponentially distributed. For n > 1 wc define 

T n := £{t>0:£ t = n} 

v n := E v [T n \T n > 0] - E v [V n \V n > O]^ 1 

so that t = J2 n >i Tn and for alH > 

V(T n >t)< V(T n > t\T n > 0) = e-*/"- 

We claim that E v [V n \V n > 0] converges to some finite limit as n — > oo. 
Then, since by Remark 8, 7r„ = a 2 X 2n ^1 — 0(X 2n ), we have v n = 

0(\~ 2n ), so that the quantities v := '}2 n>1 v n and A := — J2n>i Vn v n are 
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both finite. Define the sequence of numbers (0 n ) n >i satisfying 
e -te n /u n = I/ne (A-t)/« > n > 1 

and notice that 

= 1 



n>l n>l 

so we conclude that 
V{t > t) < v( (J {T„ > n t}) < ^ 7>(T„ > n t) < ^ e 

^n>l ' n>l n>l 

This proves the theorem for 

2 



t$ n /Vn = ve (h-t)/v 



C= \\xfve^ = N| 2 cxp{-]T v l0g vl " 11 

n>l 

^[K|K >o] 



\2n 1.2 



We check that these do not depend on a. From the definitions of v n and v, it is 
enough to show that the law of (k does not depend on a. 

The transition probabilities of £ are given by p n ,n = and p n ,m '■— for 
n 7^ to. Recall from (30) that 

(43) 7r n , m = cr 2 A 2 " k f and ^ = °" 2a2 " £ fc2 

ri=m— n 

meaning in particular that Pn.m 

does not depend on a. 
Finally, we must prove the claim. 

Consider p n ,n+r and notice that again by (43) and Remark 8, it does not depend 
on n, for n > uq. This means that, as long as > no, C behaves like a random 
walk with increment distribution 

V k 2 

L-ii:n=r i _ m 

q r :=^ — pr, ^ez. 
2-rjeJ «j 

Let pk, for fe = 0, 1, 2, . . . be a random walk on Z defined on (S, S, V) starting from 
po = Coi with increment distribution q. Since 

$>V.= £ Af = $> 2 A 2 - = 5r A"^ 

r*=— r r^— r 

and A > 1, we have that g_ r < q r whenever r > 0, so p has a positive drift. 

Now we forget for a moment the starting distribution of £ and consider only 
transition probabilities. Let H = {pk > no,Vfc > 0} and K = {(k > no,V/c > 0} 
and let n > uq. Then 

(44) V(K\( = n) = V(H\ Po = n) 

(45) V{( k ? n,Vn > l|ft = n,K) = V{p k ^n,Vn> l\p = n,H) 

Take the limit for n —> oo. Since p is a random walk with a positive drift, then (44) 
converges to 1. Hence in the limit we can drop H,K from (45) and conclude that 

lim V(Ck 7^ n,\fn > l|(o = n) = lim V{pk ^ n,Vn > l|po = n) > 

n— toe n— foo 

This in particular means that E v \y n \V n > 0] converges to some finite limit, which 
was the claim we had to prove. □ 
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The statement of Theorem 5 is about expectations, but since the decay is at 
least exponential, it can be refined to an almost sure convergence by virtue of 
Borcl-Cantelli Lemma. Proposition 16 below gives the details. 

Lemma 15. Under the same hypothesis of Theorem 5, let s > 0, then Q-a.s. 

sup\\X(t)\\ < \\X(s)\\ 

t>s 

Proof. Let x = X(s), and consider the restriction of X to the time interval [s, oo). 
Then by Theorems 4 and 2 and the ensuing remark, X is the unique strong Leray 
solution and has the property that Q- almost surely ||X(f)|| < ||x||. □ 

Proposition 16. Under the same hypothesis of Theorem 5, the total energy of the 
solution goes to zero at least exponentially fast pathwise under Q, 

lim sup — log||X(t)|| 2 < , Q-a.s. 

t^oc t (I 

Proof. Let a > 0. By Theorem 5 we can bound the probabilities 

Q(\\X(n)\\ 2 > e- an ) < Ce~- na2/f *e an , n>0 

If we take a < ct 2 //i, by Borel-Cantelli lemma there exists a r.v. M, such that 
Q-a.s. and for all n > we have \\X(n)\\ 2 < Me~ an . 

For n = 0, 1, . . . , apply Lemma 15 with s = n to get that Q-a.s. 

sup \\X(t)\\ 2 < Me~ an 

te[n,n+l) 

or 

sup \\X(t)\\ 2 e at ~ a < sup \\X(t)\\ 2 e an < M 

te[n,n+l) te[n,ri+l) 

These are countably many propositions, so Q-a.s. all of them are true, yielding 

sup||Jf(i)|| 2 e Qt < Me", Q-a.s. 
t>o 

From here the thesis follows quickly by letting a /* a 2 / ft on the rational numbers. 

□ 

To translate the almost sure statement of the above proposition to the initial 
non-linear problem, one should be able to prove the equivalence of P and Q on 
J^oo (see Remark 6). The following proposition is the key result to prove Novikov 
condition of Girsanov theorem for t — oo, which is the object of Theorem 6 below. 
A very similar statement can be found in Barbato, Flandoli and Morandin [7] and 
the proof, which is almost the same, is given here for completeness. 

Proposition 17. Under the same hypothesis of Theorem 5, let fi > be the con- 

2 

stant given there and let 9 > 0. If 9 < Tjrara > then 




Proof. Let V :— J °° ||X(£)|| 2 dt and take any v > 0. Let u > defined by 
(46) IW| 2 Q(^ >v) = Ce- U " 2 ^ 
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where p, > and C > ||a;|| 2 are the constants given by Theorem 5. Then 

poo 

vQ(V >v)< E Q (V; V>v)= / E Q (\\X{t)\\ 2 ; V > v)dt 

Jo 

POO 

< / mm(\\x\\ 2 Q(V > v);E Q (\\X(t)\\ 2 ))dt 
Jo 

< / \\x\\ 2 Q(V > v)dt + / Ce- try2 ^dt 

Jo Ju 

< u\\x\\ 2 Q(V >v)+ f ia- 2 Ce- ua2/ » 
= \\x\\ 2 Q(V > v)(u + fi<j- 2 ) 

where we used Leray property, Theorem 5 and twice equation (46). 

If Q(V > v) = for some v then V is bounded and we are done. Otherwise we 
get a lower bound on u which put into (46) gives 

1 , „2 



yielding E«(e ei/ )<oo for all 0<^p. □ 
Theorem 6. Under the same hypothesis of Theorem 5, let fi > be the constant 

\/u,\I\\\x\\ 

given there and let p :— v 2g . 2 — . If p < 1, the total energy of the solution goes to 
zero at least exponentially fast under P, both pathwise and in average, 

(47) limsupilog||X(t)|| 2 < - — , P-a.s. 

t-yoo t fl 

(48) limsup - logE p ||X(t)|| 2 < - — (1 - pf 

i->oo t fl 

Moreover P and Q are equivalent on .Too . 

Proof. Novikov condition (25) can be extended also to the case T = oo, and Propo- 
sition 17 applied with 9 = = j^z shows that it holds on (Q,Q), so P <C Q- 
The density is given by (22) with t — oo; it is a.s. positive, P and Q are equivalent. 
The first statement then follows by Proposition 16. 

To prove (48) we follow Barbato, Flandoli and Morandin [7]. Fix t > and let 
f = j^= exp{Z t - \[Z, Z] t } (see (22)). Let p,q>l with \ + \ = 1. Then 

E p \\X{t)\\ 2 = Efl(f\\X{t)f) < EQ(pf lp E9(\\X{t)f< 1 f /q 
We bound the first term by (24) and Girsanov theorem 

E«(/*)=E«(exp{pZ,-|[Z,Z] t }) 

= EQ(cxp{pZ f - \[pZ,pZ] t + PtP-pl[Z,Z} t }) 

< exp{ p(P 2 ' 1)|r| ^ l|2t }EQ(exp{ P Z t - \\pZ,pZ] t }) 

We bound the second term by Leray property and Theorem 5, 

E Q (\\X(t)\\ 2q ) < \\x\\ 2q - 2 E Q (\\X(t)\\ 2 ) < WxW^Ce'^ 1 
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Putting together the two bounds above and with some algebraic manipulations we 
get that for all p > 1 

logE p ||Y(t)|| 2 < loglMI 2 + (l - -) (p P 2 -t - ^i + log^ 



This formula can be optimized on p. The RHS attains its minimum when 

p 2 = p- 2 (l-^-lo g ( 



If t is large enough and p < 1, then this gives an acceptable value p > 1. By 
substituting this value of p and letting f->oowe get the thesis. □ 



7. Applications 

In this section we apply our general model to two important shell models of 
turbulence, namely the inviscid versions of GOY model 

(49) ^-u n = ia\ n u* n+1 u* n+2 + ibA^-iU^u*^ + icA„_ 2 <_i<_ 2 ; n > 1 
and the inviscid version of Sabra model 

(50) —u n = ia\ n u* n+1 u n+2 + ibXn-xU^^Un+i - icA„_ 2 u„_iu„_ 2 , n > 1 

In both models for n > 1, u n are complex- valued functions, A„ = A™, A > 1, a, b, 
c are real numbers with a + b + c = and we set A„ = 0, u n = for n < for 
simplicity. 

We may add multiplicative noise to both models to fall in two special cases of 
our general model (1). This must be done according to the initial requirements and 
it turns out that the proper way to add noise is for the GOY 

(51) du n = ia\ n u* n+1 u* n+2 dt + z6A Il _iw*_ 1 u* +1 dt + icA n _2U„_iW n _ 2 ^ 

+ ia\ n u* n+1 odw n - £5-A n _iti*_i °dw n -i 

and for Sabra 

(52) du n = ia\ n u* n+l u n+2 dt + ib\ n -iu* n _ l u n+1 dt - ic\ n -2U n -xU n -2dt 

+ i&i\ n u* n+1 odw n - iCTiA n _iu*-i °dw n ^i 

+ (io2A n w* +1 odw' n )* ~ ia 2 \n-iu n ~i °^n_i 

where a,ai,d- 2 are positive constants with = ^ and (u>„) n ez, {w' n ) n £Z are two 
sequences of complex- valued Brownian motions which are all independent. 

Definition 4. Given an initial condition u° £ ^ 2 (C), a Leray solution of the stochas- 
tic GOY system (51) (respectively of the stochastic Sabra system (52)) is a filtered 
probability space (fi, JFoo, {Ft)t>o, P), along with an adapted sequence (u> n )nez 
(resp. two adapted sequences (w n )„GZ and (w' n ) ne z ) of independent complex-valued 
Brownian motions, and a stochastic process u such that 

i. u — (u n ) n >i is a stochastic process on ($7, T^, (J 7 t)t>o> P) taking values 

in Z 2 (C) with continuous adapted components; 
ii. With probability 1, for all t > 0, \\u(t)\\ V 2 < \\u°\\p. 
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Hi. the following integral equation (resp. equation (54)) holds for all n > 1 
and all t > 0, 

(53) u n {t)=u° + ia\ n u* l+1 (s)u* l+2 (s)ds + i6A„_i<_ 1 (s)u* +1 (s)ds 
Jo Jo 



+ / ic\ n -2K_ 1 (s)u^ 2 (s)ds - 



t ~2 



(^n + A^_ 1 )u„(s)ds 



i — •/(.— ^ —n— i v"/ ~~n— z v~/ — / <-) 
JO JO z 

+ / iaX n u* n+1 (s)dw n (s) ~ I iCTA„_iu*_ 1 (s)dw n _i(s) 



o 



O 



(54) u n (t)=u° l + ia\ n u* l+1 (s)u n+2 (s)ds + ibX n ^ 1 u* l _ 1 (s)u n+1 (s)ds 
Jo Jo 



«cA„_2Un-l(s)M n _ 2 (s)(is - 
l] Jo 

t r t 



t X.2 , ~2 



'2a2 



(A n + \ n _ 1 )u n {s)ds 



+ idi\ n u* n+1 {s)dw n (s) - I «o'iA„_iu*_ 1 (s)<iu> n _i(s) 



o Jo 

+ / {ia 2 X n u* n+1 (s)dw' n {s))* 



ia2\n-lUn-l(s)dw' n _ 1 (s) 



Theorem 7. Given an initial condition u° £ ? 2 (C), there exists a Leray solution 
(P, u) of the stochastic GOY system which is unique in law. Moreover for allt > 0, 

p(IK*)||p < ||«°||p)>o 

and for all e > there exists t > such that 

P(|K*)||p < e) > 

Finally, if ||u°||/2 is sufficiently small, then for t — > oo, u(t) converges to zero at 
least exponentially fast both almost surely and in L? (fi; I 2 (C)) . 

Proof. All we need to do is rewrite this model in the formalism of our general 
model (1). Take d — 2, let <f) : C — > R 2 be the obvious isomorphism and let 
X n := 4>(u n ) ■= (Re(u„),Im(u„)) and x n = 4>{u°). 

The first step is to define a bilinear operator B on K 2 corresponding to (v, z) > 
iv*z* on C. For a,0,j € {1,2}, let 

(l/y/2 a + /? + 7 = 4 
B Q ' /3 ' 7 = < -1/V2 a + /? + 7 = 6 
[ otherwise 

so that it is easy to chack that for any v, z g C, <p{iv* z *) = V2B(<p(v),(f)(z)) and 
that L = L a 'P = J2 y> s B a '^ s B^' 5 is the identity. 

The second step is to choose the interactions corresponding to the GOY local 
range coupling. Since there are three terms, at least two pair of interactions are 
needed. Let I = {1,2,3,4}, r = (1 3) (2 4), 7* = {1,2} and for i e I let Bi = B 
and 



i 


Ti hi hi 


1 


1 


2 




2 


-1 


-2 


V2\- 2 c 


3 


-1 


1 


— v/2A _1 ct 


4 


1 


-1 


-y/2\- l c 



STOCHASTIC INVISCID SHELL MODELS 



21 



It is now easy to check that if we apply (ft to the sum of the first three integrals 
appearing in the RHS of equation (53), we simply get 



(55) 



^ / ki, n Bi(X n+n (s),X n+hi (s))ds 



Finally, let a := a/Va 2 + A _2 c 2 , let W — (Wj, n )j e /, ne z be a family of 2- 
dimensional Brownian motions symmetric with respect to r, let (w n )nez be a se- 
quence of independent complex-valued Brownian motions and suppose that the 
following equation holds for all n, 

(56) w n = = 2 ■ — 

Va 2 + A" 2 c 2 Va 2 + A" 2 c 2 

Then 

id-A„M* +1 (s)dw;„(s) 

= / io\ n u* n+l adWl n+2 - / iCTA„< +1 A _1 cdW^^ 
Jo Jo 



/ i'CTA n < l+1 i'adt¥ 1 2 n+2 
Jo 



+ / io\ n u* n+1 i\ 1 cdW% n _ 1 



i<Ja\ n u* n+1 (dWt n+2 -idWl n+2 )- / iaX 1 cX n u* n+1 {dWl n _ 1 -idWl n _ 1 ) 



so that we can compute (ft applied to the two stochastic integrals appearing in the 
RHS of equation (53): we obtain 

(ft\ I iaX n u* l+1 (s)dw n (s) 



o 



r-t 



= f iaV2aX n B(X n+u dW hn+2 )- f iaV2X- 1 cX n B(X n+1 ,dW 2 , n - 1 ) 
Jo Jo 

= I vh.nB l (X n+rt ,W i , n+hi ) 

and analogously 

/ ia-X n - 1 u* l _ 1 (s)dw n - 1 (s) J = V" / (jki, n Bi{X n+r ^dW^ n+ v H ) 
Jo / l= 2,3 Jo 



The last term is 



' y (A* + A 2 _i)w n (s)rfs^ = - jf y (a 2 + A~ 2 c 2 )(l + A- 2 )A 2 "X„( s )cfe 



E / Y k ^ Xn(s)ds 

iei 



We have proved that, under the assumption that (56) holds, if we apply (ft to 
equation (53), we get equation (9). But of course, given (W/ n )i e / J= i i 2, then (56) 
may be taken as a definition of w n , so existence of a Leray solution follows from 
Corollary 13. On the other hand, given w, let w be onother sequence of independent 
complex-valued Brownian motions, independent from w. Then 

fW hn+2 \ ,_ fW 3 , n+2 \ 1 / a -A-! C \ /ReK) Im( W „)\ 



W 2 , n -lJ ' \Wi, n -i) ' Va 2 + A- 2 c 2 \ A c a / V Re (*«) Im (^n) 
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defines the family W according to the requirements and to (56), so uniqueness in 
law of the Leray solution follows from Theorem 3. 

To prove the two inequalities, remember that ||X(t)|| = ||tt(£)||;2 and apply The- 
orem 5. Fix t > 0. Since E Q (\\X(t)\\) < \\x\\, then Q(||X(t)|| < ||x||) > 0, so the 
same holds for P which is equivalent to Q on J ( . 

Fix e > 0. Since E Q (\\X(t)\\) -> as t -> oo, then for t large enough Q(\\X(t)\\ > 
e) < 1, so the same holds for P which is equivalent to Q on Ft- 

Finally, to prove the last statement, we apply Theorem 6. If = ||u°||;a is 
small enough, then p < 1, so by (47) we get that P-a.s. for all e > 0, for t large 



Ht)\\l2 < e-^~ e)t 



and by (48) we get that for all e > 0, for t large 

IK*)lli W( o) = eP II«(*)IIp ^ e-(£ (1 - pf - e > □ 

Theorem 8. Given an initial condition u° € l 2 (C), there exists a Leray solution 
(P, u) of the stochastic Sabra system which is unique in law. Moreover for all t > . 

p(IK*)||p < ||«°||p)>o 

and for all e > there exists t > such that 

P(|K*)||p <e) >0 

Finally, if ||u ||p is sufficiently small, then for t — > oo, u(t) converges to zero at 
least exponentially fast both almost surely and in L (C)). 

Proof. Wc follow the same strategy as for Theorem 7, so let d, <fi, X, x, I, r, r^, 
hi and ki be defined like there. We need three new different bilinear operators on 
R 2 (below on the right) which represent the corresponding bilinear operators on C 
associated to the interactions in the Sabra model (below on the left) 

a + f3 + 7 odd 

(>. : i - n- :. B^^ 1 = B^' 1 = { -l/y/2 a = 1,0 = 1 )7 = 2 

otherwise 

a + j3 + 7 odd 

■//■; B^' 1 = \ -1/V2 a = 2,j9 = l,7 = l 

otherwise 

a + j3 + 7 odd 

('•.-) irr B2' p ' y = { -l/v/2 a = 1,^ = 2,7 = 1 

otherwise 

These Bi satisfy (3) and the corresponding Lj's are the identity. 

It is immediate to verify that if we apply 4> to the sum of the first three integrals 
appearing in the RHS of equation (54), we simply get 




(57) \] \ k l . n B i (X n+rz (s), X n+hz (s))ds 

iel Jo 



Finally, let o := b\ja = &2/ (X~ 1 c), let (iw n ) n gz an d (w' n ) ne i be two sequences 
of independent complex- valued Brownian motions and let W — (Wi, n )i£i t neZ be a 
family of 2-dimensional Brownian motions symmetric with respect to r such that 
P-a.s. W\, n = TF 3 ,„ = (j)(w n - 2 ) and W 2 , n = W 4 , n = <f>(w' n+1 ). 
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Then it also easy to verify that if we apply <f> to the sum of the four stochastic 
integrals appearing in the RHS of equation (54), we get 

(58) V / Gk l , n B l {X n+rr {s),dW l ^ +hl {s)) 

The last term is 

^("f^ (A " +tl)Un(s)(fe ) 

= - f y(a 2 + A- 2 c 2 )(l + \- 2 )\ 2n X n (s)d s = - £ f yfc^X„( S )rf S 
Jo ieI Jo 

We have proved that if we apply <f> to equation (54), we get equation (9). Then one 
concludes exactly like in Theorem 7. □ 

Remark 18. The smallness condition on ||w°||;2 can be made precise by computing 
fi as in the proof of Theorem 5. One has only to observe that in this case the 
discrete-time embedded Markov chain is a simple random walk on the positive 

A 2 — 1 

integers reflected in 1 and with positive drift A2+1 and do some computations. We 
give only the result and notice that this choice of /.i is not believed to be optimal. 
For both the stochastic GOY and Sabra models, the condition p < 1 is equivalent 
to 

\\u \\p < \/2(A - A- 1 )^ 2 - A- 2 c 2 a 2 
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